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Abstract. We observe that the classical Faulhaber's theorem on sums of odd 
^ , powers also holds for an arbitrary arithmetic progression, namely, the odd power 

^ I sums of any arithmetic progression a + b,a + 2b, ... ,a + nb is a. polynomial in na + 

n{n + 1)6/2. While this assertion can be deduced from the original Fauhalber's 
Q ■ theorem, we give an alternative formula in terms of the Bernoulli polynomials. 

^ . Moreover, by utilizing the central factorial numbers as in the approach of Knuth, 

(-H ! we derive formulas for r-fold sums of powers without resorting to the notion of 

r-refiexive functions. We also provide formulas for the r-fold alternating sums of 
powers in terms of Euler polynomials. 
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O ■ 1 Introduction 

' The classical theorem of Faulhaber states that the sums of odd powers 

B' 

^ ; 1^™"^ + 22"-^ + . . . + n^"""^ 

^ . can be expressed as a polynomial of the triangular number T„ = n{n + l)/2; See 

I Beardon [2] , Knuth [12] . Moreover, Faulhaber observed that the r-fold summation 

of n"^ is a polynomial in n{n + r) when m is positive and rri — r is even [12]. The 
classical Faulhaber theorem for odd power sums was proved by Jacobi [TT]; See 
also Edwards [1] . Let us recall the notation on the r-fold power sums: — 
n"^, and 

^'n™ = + 5^""^2'" + ■ • • + J]]""^^'"- (1-1) 

For example, ra"* = l™ + 2'" + ■ ■ ■ + ra"*, and 

n 

^2 ^"^ = ^1 1'" + ^1 2"^ + • ■ • + n"^ = + 1 - 

i=l 

For even powers, it has been shown that the sum 1^"* + 2^™" + ■ ■ • + ra^*" is a poly- 
nomial in the triangular number T„ multiplied by a linear factor in n. Gessel and 



1 



Viennot [H] had a remarkable discovery that the alternating sum 1 
is also a polynomial in the triangular number T„. 
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Faulhaber's theorem has drawn much attention from various points of view. 
Grosset and Veselov [7] investigated a generalization of the Faulhaber polynomials 
related to elliptic curves. Warnaar [15], Schlosser [I^ and Zhao and Feng [16] 
studied the g-analogues of the formulas for the first few power sums. Garrett 
[5] found a combinatorial proof of the formula for sums of g-cubes. Guo and 
Zeng [S] obtained the g-analogue formula in the general case. Furthermore, Guo, 
Rubey and Zeng |9] have shown that the g-Faulhaber and g-Salie coefficients are 
nonnegative and symmetric in a combinatorial setting of nonintersecting lattice 
paths. 

In this paper, we first formulate Faulhaber's theorem in a more general frame- 
work, that is, in terms of power sums of an arithmetic progression. Given an 
arithmetic progression: 

a + 6, a + 26, . . . , a + nh, 

Faulhaber's theorem implies that odd power sums of the above series are poly- 
nomials in na + n{n + 1)6/2. In particular, an odd power sum of the first n odd 
numbers 

^2m-l ^ 32m-l + . . . + - l)^™-^ 

is a polynomial in n^, and the sum 

is a polynomial in the pentagonal number n(3n — l)/2. 

Because of the relation (a + hi)™' = b''^{a/b + z)"*, there is no loss of generality 
to consider the series 

x + l,x + 2,...,x + n. (1.2) 

Let 

X = n{n + 2x + l) (1.3) 
be the sum of the sequence x + l,x + 2, . . . ,x + n. The the power sums 

S2m-i = ix + l)2™-i + (x + 2)2™-i + . . . + (x + n)2— 1 

is a polynomial in A. For exmaple, 

1 1 
^5 = — + —(6x^ + 6x -1)A2 + -(3x^ + 6x^ + 2x2 -x)A; 
6 12 6 

A^ 1 1 
Sj = _ + -(3x2 + 3x-l)A^ + — (9x^ + 18x^ + 3x2 -6x + 1)A2 
8 6 12 

+ -(3x^ + 9x^ + 6x^ - 3x^ - 2x^ + x)A. 
6 
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It should be noticed that the above more general setting of Faulhaber's theo- 
rem can be deduced on the the original version of Faulhaber's theorem. When x 
is a positive integer, we have the relation 

n n+x X 

1=1 i=l i=l 

By Faulhaber's theorem, Ym=i ^^"^ YTi=i polynomials in (n+x)(n+x+l) 

and x{x + 1), respectively. Using the following simple but important relation 

(n + x){n + X + 1) = n{n + 2x + 1) + x{x + 1), (1.4) 

we see that 

[(^ + + a; + 1)]^ _ + 1)]^ = ^ Q [n{n + 2x + l)f[x{x + (1.5) 

which is a polynomial in n{n + 2x + 1). Clearly, one sees that the above assertion 
holds for all real numbers x. 

Although in principle Faulhaber's theorem is valid for any arithmetic pro- 
gression, from a computational point of view it still seems worthwhile to find 
a formula for the coefficients in terms Bernoulli polynomials. The main result 
of this paper is an approach to the study of the r-fold sums of powers without 
resorting to the properties of r-reflective functions as in the approach of Knuth 
|12j . In the last section, we obtain formulas for the r-fold alternating sums of 
powers in terms of the Euler polynomials. 



2 An Alternative Formula 

In this section, we give an explicit formula for the coefficients regarding Faul- 
haber's theorem for the series x + 1, x + 2, . . . , x + n, which reduces to an al- 
ternative formula to the Gessel-Viennot formula when setting a; = 0. We first 
recall some basic facts about Bernoulli polynomials Bn{x) which are defined by 
the following generating function: 

^ g„(x)r _ te^' 

^ n\ "e*-l- ^ ' 

The power sums of the first n positive integers can be expressed in terms of Bi{x): 



= -^(B^+iin + 1) - S™+i(l)). 

1=1 

Moreover, we have 

n 

Y^{x + if^-^ = — {B2m{x + n + 1) - B2m{x + 1)) . (2.2) 



2m 

i=l 
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The Bernoulli numbers -B„ are given by -B„ = -B„(0). Note that the Bernoulli 
polynomials satisfy the following relations 

Br,{x + 1) - Br^ix) = nx'^-\ (2.3) 
Bn{l-x) = (-irs„(x), (2.4) 

^B^{x) = n5,_i(x), (2.5) 
B^{x + y) = X;(^')5,(x)y'^-\ (2.6) 

The evaluation of Bernoulli polynomials at 1/2 is of special interest. For n > 0, 
we have 

E2„+i(l/2) = 0, 52„(l/2) = (2^-2" - l)B2n. (2.7) 
From (12. 6p and (12.71) . we can deduce the following form of Faulhaber's theorem. 

Theorem 2.1 Let A = n{n + 2x + 1). Then we have 

n m 

Y,{^ + rf-^ = Y,Ft\x)X\ (2.8) 



i(m) 

i=l k=l 

where 



Proof. From the binomial expansion (12. 6p . we get 

52™(a; + n+l) = ^ (^^^^52^_, (^0 • + n + . (2.10) 
It follows from and ([22D that 

- ")^"-' - L t ) (5) ■ ^ " 5) " - ^1)- 

Since 

. + „ . 1)^' . (a + (. + . t Q (. + (2.11) 

we immediately get (12. 9p for k ^ 1. For = 0, we have 

This completes the proof. I 
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The above formula can be viewed as an alternative form of the formula of 
Gessel and Viennot [B] for the coefficients A^j^^ = F^"^\o): 

^—^ \m — k — J J \ J J m — k + J 

Note that Bq = 1 and Bi = —1/2 are used in the above formula whereas 
the formula (12.91) does not involve Bi. The equivalence between the formulas 
for ^^^""^(0) and can be established via the following generating function for 
the coefficients fI"^~^^\x). The proof is analogous to that given by Gessel and 
Viennot [6] for the case a; = 0. 

Theorem 2.2 We have 



oo oo „,2m+i coshy./ (x + i)2 + t - coshy(x + i) 

l^l^^k Wf(2m + 1)!" 2sinh(f) ■ ^ ' 

m=0 fc=l \ < J \2' 

Similarly, the theorem of Gessel and Viennot on the alternating sums of even 
powers of the first n natural numbers can be extended to an arithmetic progression 
x + l,x + 2,...,a; + n. It turns out that the Euler polynomials play the same role 
as the Bernoulli polynomials for sums of odd powers. 

The Euler polynomials En{x) are defined by 

oo in Op^t 

n=0 

The following expansion formula holds: 

En{x + y) = J2(f)Ek{x)y--'. (2.14) 

fc=0 ^ ^ 

For positive even number n, we have -E„(l) =0 . The Euler numbers E„, and the 
Euler polynomials are related by 

E2n+i = 0, K = 2"E„(l/2), n>0. (2.15) 
Theorem 2.3 Let A = n{n + 2x + 1). Then we have 

n m 

J](-l)-(x + z)^- = J]Gf)(x)A^ 

i=l k=0 



where 



lv^/2m\/A_ /1\ / 1^'~^ 



i=k 



Gt'\x) = ]^{l-{-ir)EUx + l). 
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The generating function for G^IJ^\x) is given below, which is a straightforward 
extension of the formula of Gessel and Viennot [6]. 

Theorem 2.4 We have 



^ , ^ 7,2m cosh (X + + t 

^2m)! 2cosh(|) 



3 r-Fold Sums of Powers 

In this section, we derive a formula for the r-fold sums of powers of the series 
X + l,x + 2, . . . ,x + n in terms of the central factorial numbers as used in the 
approach of Knuth [12j. A key step in our approach is the r-fold summation 
formula for the lower factorials. It can be seen from our formula that if r and m 
have the same parity then the r-fold power sum Yl^i''^ + x)'" is a polynomial in 
n{n + 2x + r) plus a term that vanishes when a; = 0. The cases when r and m 
have different parities can be dealt with some care, and the details are omitted. 

Recall the notation for the lower factorials {x)k = x{x — l)---{x — k+l). The 
definition of the central factorials 213], is given by 

= x{x + k/2 - l)k-i. 

The central factorial numbers T{m, k) are determined by the following relation: 

m 

a;™ = ^T(m,A;)x['=l, m > 1. (3.1) 
fc=i 

Note that T(m, k) = when m — k is odd. In particular, we need the following 
relation 

m 

^2m-i ^ Y^T{2m,2k){x + k-l)2k-i- (3.2) 

k=l 

We first give a formula for the r-fold sums of lower factorials. From the 
recursive definition fll.ip of r-fold summations, we have 

Y^\x + n)i= + (3-3) 

The above multiple summation can be simplified to a single sum. 
Theorem 3.1 We have 

{x + n + r)i^r sr^ fn + r- i-l\{x + 



(/ + r)r V ^ ~ ^ 7 + ^) 
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Proof. Setting 

= y^(x + n + 1)^+1 
gives F{i) — F{i — 1) = (x + Hence we get 



1^*1^12 41=1 



^ 'x + i2 + l)/+i-7^(a; + l)«+i. (3.4) 



/ + r l + l 

Iterating (13. 4p r — 1 times and using the following identity 



E 

i=l 

we obtain the desired formula. 



From Theorem 13.11 and the relation (13.21) . we derive two formulas for the r-fold 
sums of the m-th powers when r and m have the same parity. 

Theorem 3.2 Form > 1, we have 
E"^\x + n)-- = ET(2m,2fc)^(^ + ^ + ^ + 2r),.^^ 



{2k + 2r)2.+i 

2k+i-l 



k=l 

n + 2r — A (.x + A; + z — 1 



^^\2r-%\\) (2A; + z-l) 

We remark that the second summation in the above formula vanishes when 
X = 0, and the lower factorial (n + + 2r)2A;+2r can be rewritten as 

]^(n + A; + 2r + 1 - i)(n - A; + i) = ]^[n(n + 2r + 1) - (A; + 2r - i + 1)(A; - i)]. 

i=l i=\ 

Hence we obtain Faulhaber's theorem for the (2r + l)-fold sums of odd powers of 
the first n positive integers. 

Applying Theorem 13.11 together with the following relation 
(x + n)(x + n + A; - l)2fc-i = ^(x + n + A;)2fc + ^(x + n + A; - 1)2^, (3.5) 
we arrive at the following formula for the (2r)-fold summation of even powers. 

Theorem 3.3 Form > 1, we have, 

y-\x + nf- = V T(2m, 2A:) | + ^ + + " + ^ + ' ^^^^-^""^ 
i^^i I (2«^ + 2r)2r 

2r 



El n + 2r - i - 1\ {2x + i){x + k + i - l)2fc+i-i 
V 2r-i / 



2(2A; + z) 
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Setting a; = in the above formula, the second summation vanishes. One 
sees that the (2r)-fold summation becomes a polynomial in n{n + 2r) because 
{n + r){n + k + 2r — l)2k+2r-i can be rewritten as 



k+r k+r 



Y[{n + k + 2r - i){n - k + i) = JJ[n(n + 2r) + {k + 2r - k)]. 



i=i 1=1 



4 r-Fold Alternating Sums of Powers 

In this section, we investigate the r-fold alternating sums of powers. Following 
the notation of Faulhaber, we define 

j2\-mn+xr:= Yl (-ir(^i+^r- (4-i) 

l<ii<---<ir<n 

We will show that the 2r-foId alternating sum of even powers ^^'^( — is a 
polynomial in n{n + r). For other cases concerning the parities of r, m, we will 
outline the results without proofs. 

Define 

E f- .)E,,{x,)---Eaxr). (4.2) 

The following lemma holds. The proof is based on induction and is omitted. 
Lemma 4.1 Let r , m he positive integers. Then 

= (-l)"2-'i?M (^1/2, . . . , 1/2, (a; + n + r/2 + 1/2)) 

+ E ^ ' : r ') 2"'^^'^ (V2> ■ ■ ■ ' 1/2, X + + l)/2) . 

k=l ^ ^ 



(k) 

We now give recursive formulas for i^g^ in order to compute the multiple sums 
in the above lemma. 

Lemma 4.2 Let k, m he positive integers. Then e'^^ [1/2, . . . , 1/2, [k + l)/2] 
equals 

E(2i)EQ(-l)'^'^»'(V2....,l/2), (4.3) 

and (1/2, ... , 1/2, {k + l)/2) equals 

E (si + l) E ( •) (-1)'+'^^^^:." (1/2, .... 1/2, 1) . (4.4) 
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Proof. From the generating function of Em{x), one sees that the exponential 
generating function of Em'^ ^1/2, . . • , 1/2, {k + equals 

E^«(l/2,....lA(*+l)/2)i^^ ^ . 

m=0 ^ 

Observe that 




Denote the exponential generating functions of and En{l/2) by A{t) and 

B{t), respectively. Then the above relation implies that = 2A{t) — B{tY, 

which yields 

A{t) = 1 - ^l-B{tf. 
Taking the fc-th power, we obtain 

Since i?2m+i(l/2) = and i?2m = for m > 1, equating the coefficients of 
both sides of above identity yields the desired recurrence relations. I 

The following theorem is concerned with 2r-fold alternating sums of even 
powers. 

Theorem 4.3 Let r, m be positive integers. Then the 2r-fold alternating sum 
IS of the form (-l)"F(n(n + 2r)) + G(n(n + 2r)), where F and G 
are polynomials of degree m and r — 1 respectively. 

Proof. Applying Lemma WTLl with 2r replaced by 2, m replaced by 2m, and setting 

X = 0, we get 

= (-l)"2-2'-i?£) (l/2, . . . , 1/2, (n + r + 1/2)) (4.5) 

+ E ^ 2r I r 2"^^- • • • ' + • 

fc=i ^ ^ 

In the expansion of 04.51) indexed by 2ii + 2i2 + - ■ ■ + 2z2r = 2m, each term contains 
a factor of the form E2i^{n + r + 1/2). According to (12.141) . we find that 
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Since (n + r)^'^ = {n{n + 2r) + r^)'^, it is easily seen that (14.51) is a polynomial 
in n{n + 2r) of degree m. We need to show that fl4.6p is also a polynomial in 
n{n + 2r). Applying (14.31) . we find that (14.61) equals 



fc=l ^ / j=0 ^ ^ j=0 ^-^^ 

r 

= E(-l)'^2^^(V2,...,l/2) 



X 



Using the identity 

n/2 



we deduce that the second sum in (14.71) equals 



k—j 



In view of the following relation, 

E 



n — k\ f m + k\ fm + n + 1 
3 



the above sum (14.91) simplifies to 



n + 2r — j — l\n + r 



-2j 



(4.10) 



2r — 2j — 1 J r — j 
Now, {n + r){n + 2r — j — l)2r-2j-i can be rewritten as 

r-j r-j 

Y[{n + i+j){n + 2r - i - j) = JJ[n(n + 2r) + (2r - j - + j)]. 

i=l 1=1 

which is a polynomial in n{n + 2r) of degree r — j. Since e!^{-) = 0, (14. 6 p is a 
polynomial in n{n + 2r) of degree r — 1. This completes the proof. I 

For the remaining three cases with respect to the parities of r and m, we have 
the following theorem. The proof is omitted. 
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Theorem 4.4 For r,m >0, we have 

J2 (-1)"^^"^ (4.11) 
= (-l)"F«(n(n + 2r + 1)) + (2n + 2r + l)Gi'i(n(n + 2r + 1)), 

^ (-1)V"+' (4.12) 
= ( — 1) (n + r)Fj^\n{n + 2r)) + {n + r)G'fi(n(n + 2r)), 

J2 i-lTn^""^^ (4.13) 
= (-l)"(n + r)Fj^\n{n + 2r + 1)) + G^^\n{n + 2r + 1)), 

where Fm{x) and Gr\x) {i — 1, 2, 3) stand for polynomials of degrees m and r 
respectively, and G^^i — for i — 1,2. 
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